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We measured the angular rotation and proper motion of the Triangulum Galaxy 
(M33) with the Very Long Baseline Array by observing two HoO masers on 
opposite sides of the galaxy. By comparing the angular rotation rate with the 
inclination and rotation speed, we obtained a distance of 730 ± 168 kilopar- 
secs. This distance is consistent with the most recent Cepheid distance mea¬ 
surement. M33 is moving with a velocity of 190 ± 59 km s _1 relative to the 
Milky Way. These measurements promise a new method to determine dynam¬ 
ical models for the Local Group and the mass and dark matter halos of M31, 

M33 and the Milky Way. 

Measuring the proper motion and geometric distances of nearby galaxies has been a long¬ 
standing problem. As part of a famous debate about the nature of galaxies, van Maanen - an 
experienced observer — claimed in 1923 to have measured a large proper motion and angular 


1 


rotation rate for the Triangulum Galaxy (M33) on photographic plates separated by 12 years (0. 
This was proven wrong by Hubble through the discovery of Cepheids in M33 that showed a 
large distance (0. This pushed the detection of galaxy proper motions beyond the capabilities 
of past experiments. Yet, galaxy proper motions are important for many astrophysical issues, of 
which two are addressed in this report. 

First, measuring accurate distances is of great importance to all fields of astrophysics, from 
stellar astronomy to cosmology. The calibration of most standard candles used for measur¬ 
ing extragalactic distances is tied directly or indirectly to the distance to one galaxy, the Large 
Magellanic Cloud (LMC), which remains controversial ©t?}). Hence, it is important to ob¬ 
tain geometric distances to nearby galaxies in which well understood standard candles can be 
studied. This allows independent calibration and verification of the extragalactic distance scale. 

Another important issue is the distribution of luminous and dark matter in the local universe. 
The problem when trying to derive the gravitational potential of the Local Group of galaxies © 
is that usually only radial velocities are known from the Doppler effect and statistical approaches 
have to be used (00. The proper motions of some nearby galaxies in the Milky Way subgroup 
have been obtained from comparing historic photographic plates (00, but a confirmation of 
these measurements will require decades. With Very Long Baseline Interferometry (VLBI) 
and phase-referencing techniques, the proper motions for galaxies within the Local Group, on 
scales of tens of microarcseconds per year, can now be measured. The most suitable strong 
and compact radio sources in the Local Group for such a VLBI experiment are the strong H 2 0 
masers in M33 and in the galaxy IC 10 mm 

We observed two HII regions (f72l ) in M33 with known H 2 0 maser activity - M33/19 and 
IC 133 - with the NRAO Very Long Baseline Array (VLBA) eight times between March 2001 
and January 2004. M33/19 is located in the south eastern part of M33, whereas IC 133 is located 
in the north east of M33 ( Fig.l ). Our observations are grouped into four epochs, each compris- 
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Table 1: Details of the observations: Observing date, observation length t obs , beam size 6 and 
position angle PA. 


Epoch 

Date 

tobs 

9 [mas] 

PA[°] 

I 

2001/03/27 

10 

0.88x0.41 

164 

I 

2001/04/05 

10 

0.86x0.39 

169 

II 

2002/01/28 

10 

0.62x0.33 

176 

II 

2002/02/03 

10 

0.71x0.33 

175 

III 

2002/10/30 

10 

0.87x0.38 

171 

III 

2002/11/12 

10 

0.84x0.36 

165 

IV 

2003/12/14 

12 

0.85x0.36 

159 

IV 

2004/01/08 

12 

1.15x0.47 

164 


ing two closely spaced observations to enable assessment of overall accuracy and systematic 
errors (Table Q}. The separations of the two observations within each epoch were large enough 
that the weather conditions were uncorrelated, but small enough that changes in the positions 
were negligible during this time. 

We observed four 8 MHz bands, in dual circular polarization. The 128 spectral channels 
in each band yielded a channel spacing of 62.5 kHz, equivalent to 0.84 km s -1 , and covered a 
velocity range of 107 km s -1 . The observations involved rapid switching between the phase- 
calibrator J0137+312, which is a compact background source with continuum emission, and 
the target sources IC 133 and M33/19 in the sequence J0137+312 - IC 133 - J0137+312 - 
M33/19 - J0137+312. With source changes every 30 seconds, an integration time of 22 seconds 
was achieved. The background source was unresolved in all epochs and was assumed to be 
stationary on the sky. Because J0137+312 was separated by only 1° on the sky from the masers, 
we obtained a precise angular separation measurement for all sources. 

The data were edited and calibrated with standard techniques in the Astronomical Image 
Processing System (AIPS) as well as zenith delay corrections (HU) . The masers in IC 133 and 
M33/19 were imaged with standard techniques in AIPS. In IC 133, we detected 29 distinct 


3 








emission features in position and the spatial distribution was very similar to earlier observations 
mmm All components were unresolved. In M33/19, we detected eight maser features. 
Two features were separated by less than a beam size and blended together. These two features 
were fit by two elliptical Gaussian components simultaneously. All other features were fit by a 
single elliptical Gaussian component. 

The maser emission in M33/19 and IC 133 is variable on time scales less than 1 year. Be¬ 
tween the epochs, new maser features appeared and others disappeared. However, we were able 
to detect and follow the motions of four features in M33/19 and six features in IC 133 over all 
four epochs. The feature identification was based on the positions and radial velocities of the 
maser emission. Each feature was usually detected in several frequency channels. A rectilinear 
motion was fit to each maser feature in each velocity channel separately. We discarded fits with 
reduced y 2 larger than 3 as they are likely affected by blending or component misidentification. 
All features showed consistent motions within their errors (2a). The accuracy and number of 
measured motions was not adequate to model the internal dynamics of the IC133 and M33/19 
regions (such as outflow) as was done in earlier observations (17511701) . 

We then calculated the variance weighted average of all motions. This yielded an average 
motion of the maser features in M33/19 of 35.5 ± 2.7 /ras year -1 in right ascension and —12.5 
± 6.3 /ras year -1 in declination relative to the background source. For IC 133 we get an average 
motion of 4.7 ± 3.2 /ras year -1 in right ascension and —14.1 ± 6.4 /ras year -1 in declination. 

We also calculated the average position of all maser features for each observation (Fig. 2 
and 3). We used the individual fits for each maser feature to remove a constant position offset 
for each maser feature. We used the position offsets at the time 2002.627, which is in the middle 
of our observations. We then calculated the variance weighted average of the positions of all 
detected features. A fit of a rectilinear motion to these average positions yielded motions of 
37 ± 5 /ras year -1 in right ascension and —13 ± 6 /ras year -1 in declination for M33/19. For 
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IC 133 we obtained a motion of 3 ± 3 //as year -1 in right ascension and —13 ± 10 //as year -1 in 
declination. This is consistent with the variance weighted average of all maser feature motions 
and suggests that the systematic internal motions within the two regions (such as outflow) are 
probably not a substantial source of bias. The difference in the maser position of the two 
closely spaced observations within each epoch was used to estimate the accuracy of the position 
measurements. The average position error over the four epochs was 7.7 //as in right ascension 
and 9.8 //as in declination. 

The relative motions between M33/19 and IC 133 are independent of the proper motion of 
M33 and any contribution from the rotation of the Milky Way. Knowing the rotation curve and 
inclination of the galactic disk we can predict the relative angular motion of the two masing 
regions as a function of distance. The rotation of the HI gas in M33 has been measured and the 
measured velocities were fit with a tilted-ring model CH- We used this model of the rotation 
of M33 to calculate the expected transverse velocities of M33/19 and IC 133. For M33/19, we 
expect a motion of 42.4 km s -1 in right ascension and —39.6 km s -1 in declination. For IC 133 
we expect —64.0 km s -1 in right ascension and —74.6 km s -1 in declination. This gives a 
relative motion of 106.4 km s -1 in right ascension and 35 km s -1 in declination between the 
two regions of maser activity. 

The radial velocities of the H 2 0 masers in M33/19 and IC 133 and the HI gas at the same 
positions agree (< 10 km s -1 ). This suggests that the maser sources are moving with the HI 
gas in the galaxy. However, although the rotation model and the radial velocity of the HI gas at 
the position of IC 133 is also consistent (< 5 km s -1 ), there is a difference of ~ 15 km s -1 at 
the position of M33/19. This indicates the presence of another tilt in the disk that is not covered 
by the model. Because of this uncertainty in the rotation model, we conservatively assume a 
systematic error of 20 km s -1 in each velocity component for the relative velocity of the two 
maser components. 
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Comparing the measured angular motion of 30.8 ± 4 //as year 1 in right ascension with the 
expected linear motion of 106 ± 20 km s _1 , one gets a geometric distance of 


D = 730 ± 100 ± 135 kpc, 

where the first error indicates the statistical error from the proper motion measurements while 
the second error is the systematic error from the rotation model. After less than three years 
of observations, the uncertainty in the distance estimate is dominated by the uncertainty of the 
rotation model of M33. 

Within the errors the geometric distance of 730 ± 100 ± 135 kpc is consistent with recent 
Cepheid and tip of the red giant branch (TRGB) distances of 802±51 kpc and 794±23 kpc 
respectively (/T: 73). It also agrees with a geometric distance estimate of 800±180 kpc UTk . 

The observed proper motion v prop of a maser (e.g., M33/19 or IC 133) in M33 can be de¬ 
composed into three components v prop = v rot + v Q + u M33 . Here v rot is the motion of the maser 
due to the internal galactic rotation in M33 and v Q is the apparent motion of M33 caused by the 
rotation of the Sun around the Galactic center. The last contribution, vm 33 , is the proper motion 
of M3 3 relative to the Milky Way. 

The motion of the Sun can be decomposed into a circular motion of the local standard of rest 
(LSR) and the peculiar motion of the Sun. The peculiar motion of the Sun has been determined 
from Hipparcos data (tTffl) to be U 0 =10.00±0.36 km s -1 (radially inwards), V 0 =5.25±0.62 
km s -1 (in the direction of Galactic rotation) and W 0 =7.17±0.38 km s -1 (vertically upwards). 
The IAU adopted LSR moves with a velocity of 220 km s -1 towards a Galactic longitude of 
l = 90° and latitude of b = 0° (12ffb . New VLBI measurements of the proper motion of Sgr A*, 
the compact radio source at the Galactic center, indicate a slightly higher circular velocity of the 
LSR of 236±15 km s -1 , for a distance of the Sun from the Galactic center R 0 =8 kpc, where the 
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uncertainty is dominated by the uncertainty in the distance to the Galactic center (TDT) . Using a 
LSR velocity of 236±15 km s -1 plus the peculiar velocity of the Sun from d, the motion of 
the Sun causes an apparent proper motion of a & = 52.5 ± 3.3 /ras year -1 in right ascension and 
S Q = —37.7 ± 2.4 /ras year -1 in declination, assuming a distance of 730 kpc and the Galactic 
coordinates of M33 (l = 133.6°, b = —31.3°). 

Using the rotation model of m, the contribution from the rotation of M33 (for IC 133) is 
a rot = —18.5 ± 6 /ras year -1 in right ascension and 8 rot = —21.6 ± 6 fias year -1 in declination. 
Here, we assumed again an uncertainty of 20 km s -1 for the rotation velocity and a distance of 
730 kpc. Combining these velocity vectors, we get the proper motion of M33: 


^M33 &prop &rot 

= ( 4.7 ± 3.2 + 18.5 ±6 — 52.5 ± 3 . 3 )^ 

— 29.3 ± 7 . 6 -^ = —101 ± 35 — 

year s 

and 

8m 33 8p r0 p 8 ro t 5 q 

= (— 14.1 ± 6.4 + 21.6 ± 6 + 37.7 ± 2 . 4 )-^ 

45.2 ± 9 . 1 -^ = 156 ± 47 —. 

year s 

The transverse velocity changes by less than 5 km s -1 if we use the TRGB distance of 
794±23 kpc <0 for this analysis. Finally, the systemic radial velocity of M33 is —179 km s 1 
d. The radial component of the rotation of the Milky Way towards M3 3 is —140 ± 9 km s 1 . 
Hence, M33 is moving with —39 ± 9 km s -1 towards the Milky Way. This gives now the three 
dimensional velocity vector of M33 ( Fig.4 ). The total velocity of M33 relative to the Milky 
Way is 190 ± 59 km s -1 . 

For Andromeda (M31), only one component of the three-dimensional velocity vectoris 
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known, the radial velocity of 116 km s _1 (301 km s -1 systemic velocity minus 185 km s -1 
contribution from solar motion) towards the Milky Way. However, the Milky Way is possibly 
falling towards M31, because there are no other large galaxies in the Local Group to generate 
angular momentum through tidal torques (I27T) . Following this argument, we assume a proper 
motion of 0 for M31. The geometry of the Andromeda subgroup depends on the relative dis¬ 
tance between M31 and M33. Thus, it is crucial to use distances for the two galaxies that have 
similar systematic errors, and we used the TRGB distances to M33 and M31 (li#11221) . Then the 
angle between the velocity vector of M3 3 relative to M31 and the vector pointing from M3 3 
towards M31 was 30°±15°. For an angle of 30°, only elliptical orbits with eccentricities of 
e>0.88 are allowed. For the largest allowed angle of 45°, the eccentricities are restricted to 
e>0.7. This eccentricity limit weakens if the proper motion of M31 is non-negligible and, for a 
motion of >150 km s -1 , any eccentricity is allowed. 

If M33 is bound to M31, then the relative velocity of the two galaxies must be smaller than 
the escape velocity. This gives - for a zero proper motion of M31 - a lower limit for the mass of 
M31 of 1.2 x 10 12 Mg. A substantial proper motion of M31 could reduce or increase the relative 
velocity and the lower mass limit of M31 . On the other hand, the dynamical friction of M31 
on M33 indicates that M31 cannot have a very massive halo of more than ~ 10 12 M 0 unless 
the orbit of M3 3 has a low eccentricity. Otherwise, the dynamical friction would have led to 
a decay of the orbit of M33 This agrees with a recent estimates of 12.3^g 8 x 10 11 M 0 
derived from the three-dimensional positions and radial velocities of its satellite galaxies (C??lI . 

More than 80 years after van Maanen’s observation, we have measured the rotation and 
proper motion of M33. These measurements provide a method to determine dynamical models 
for the Local Group and the mass and dark matter halo of Andromeda and the Milky Way. 
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Figure 1: The positions of two regions of maser activity in M3 3. Predicted motions due to 
rotation of the HI disk are also shown. Image courtesy of Travis Rector (NRAO/AUI/NSF and 
NOAO/AURA/NSF), David Thilker (NRAO/AUI/NSF) and Robert Braun (ASTRON). 
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Figure 2: Average position of the maser M3 3/19 in (top) right ascension and (bottom) declina 
tion relative to a background source. 
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Figure 3: Average position of the masers in (top) IC 133 in right ascension and (bottom) decli¬ 
nation relative to a background source. 
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Figure 4: Schematic view of the Local Group with the space velocity of M33 and the radial 
velocity of M31. The blue cross marks the position of the Local Group Barycenter (LG BC) 

<EH>. 
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